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ABSTRACT 

The theoretical basis and architecture of the conceptual design tool NDARC (NASA Design and 
Analysis of Rotorcraft) are described. The principal tasks of NDARC are to design (or size) a 
rotorcraft to satisfy specified design conditions and missions, and then analyze the performance of the 
aircraft for a set of off-design missions and point operating conditions. The aircraft consists of a set of 
components, including fuselage, rotors, wings, tails, and propulsion. For each component, attributes 
such as performance, drag, and weight can be calculated. The aircraft attributes are obtained from the 
sum of the component attributes. NDARC provides a capability to model general rotorcraft 
configurations, and estimate the performance and attributes of advanced rotor concepts. The software 
has been implemented with low-fidelity models, typical of the conceptual design environment. 
Incorporation of higher- fidelity models will be possible, as the architecture of the code accommodates 
configuration flexibility, a hierarchy of models, and ultimately multidisciplinary design, analysis and 
optimization. 


INTRODUCTION 

The objectives of rotorcraft design work in a government 
laboratory are to support research and to support rotorcraft 
acquisition. Research activities require a robust design 
capability to aid in technology impact assessments and to 
provide system level context for research. At the applied 
research level, it is necessary to show how technology will 
impact future systems, and justify the levels of investment 
required to mature that technology to an engineering 
development stage. Design provides one avenue to 
accomplishing these objectives. The Department of 
Defense (DoD) acquisition phases requiring rotorcraft 
design work include concept exploration, concept 
decision, concept refinement, and technology 
development. During these acquisition phases, performing 
quantitative evaluation and independent synthesis of a 
wide array of aircraft designs is typically necessary, in 
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order to provide the foundation for specification and 
requirement development. 

Rotorcraft conceptual design consists of analysis, 
synthesis, and optimization to find the best aircraft 
meeting the required capabilities and performance. A 
conceptual design tool is used for synthesis and analysis of 
rotorcraft. These tools historically have been low fidelity 
for rapid application. Such sizing codes are built around 
the use of momentum theory for rotors, classical finite 
wing theory, a referred parameter engine model, and semi- 
empirical weight estimation techniques. The successful 
use of a low-fidelity tool requires careful consideration of 
model input parameters and judicious comparison with 
existing aircraft to avoid unjustified extrapolation of 
results. 

The helicopter industry has proprietary conceptual design 
tools, including PRESTO (Bell Helicopter), RDM 
(Sikorsky Aircraft), and HESCOMP and VASCOMP 
(Boeing). Until now the tools available to the U.S. 
government have been characterized by out-of-date 
software and limited capabilities. Examples are 
HESCOMP and VASCOMP (the versions developed by 
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Boeing in the 1970s), and RC (developed by the U.S. 
Army AFDD in the 1990s). 

NASA, with support from the U.S. Army, conducted in 
2005 the NASA Heavy Lift Rotorcraft Systems 
Investigation (ref. 1), focused on the design and in-depth 
analysis of rotorcraft configurations that could satisfy the 
Vehicle Systems Program (VSP) technology goals. The 
VSP technology goals and mission were intended to 
identify enabling technology for civil application of heavy 
lift rotorcraft. The goals emphasized efficient cruise and 
hover, efficient structure, and low noise. The requirements 
included carrying 120 passengers over a 1200 nm range, 
350 knots at 30,000 ft altitude. The configurations 
considered included the Large Civil Tiltrotor (LCTR), 
Large Civil Tandem Compound (LCTC), and Large 
Advancing Blade Concept (LABC). This project is an 
example of the role of a rotorcraft sizing code within a 
government laboratory. The design tool used was the 
AFDD RC code. The project illustrated the difficulties 
adapting or modifying a legacy code for configurations 
other than conventional helicopters and tiltrotors. 

Since 2005, there have been numerous other joint 
NASA/U.S. Army investigations of advanced rotorcraft 
concepts, covering conventional tiltrotors and helicopters, 
slowed-rotor compound helicopters (ref. 2), a tilting- 
tandem concept, heavy-lift slowed-rotor tiltrotors (ref. 3), 
lift-offset rotor concepts (ref. 4), and a second generation 
large civil tiltrotor (LCTR2, ref. 5). These design projects 
have gone well beyond the conventional boundaries of the 
conceptual design process, combining high-fidelity 
analyses (including rotorcraft comprehensive analysis, 
computational fluid dynamics, and structural analysis) 
with the conceptual design tool. This approach has been 
required because of the increasing sophistication of the 
requirements and the technology, and the increased level 
of certainty needed to differentiate between system 
concepts. 

Based on this experience, a new conceptual design tool 
has been developed to support future needs of the NASA 
Subsonic Rotary Wing project and the U.S. Army AFDD 
Advanced Design Office: NASA Design and Analysis of 
Rotorcraft (NDARC). The software development started 
in January 2007, and the initial code release occurred in 
May 2009. This paper summarizes the NDARC theoretical 
basis and architecture; the complete description is in 
reference 6. A companion paper (ref. 7) presents 
validation and demonstration results from the NDARC 
development. 


REQUIREMENTS 

Based on the recent experience of NASA and AFDD at 
Ames Research Center, the following requirements were 
defined for the new conceptual design tool. 

The principal tasks of the tool are to design rotorcraft to 
meet specified requirements, and then analyze the 
performance of the aircraft for a set of flight conditions 
and missions. Multiple and flexible design requirements, 
from specific flight conditions and various missions, must 
be accommodated in the sizing task. The aircraft 
performance analysis must cover the entire spectrum of 
the aircraft capabilities, and component performance and 
engine models must cover all operating conditions. A 
general and flexible capability to define conditions and 
missions is required. 

For government applications and to support research, it is 
important to have the capability to rapidly model general 
rotorcraft configurations, including estimates of the 
performance and attributes of advanced rotor concepts and 
the capability to model the impact of new technology at 
the system and component level. In such an environment, 
software extensions and modifications will be routinely 
required to meet the unique requirements of individual 
projects. Complete and thorough documentation of the 
theory and its software implementation is essential, to 
support development and maintenance and to enable 
effective use and modification. 

The code architecture must accommodate configuration 
flexibility and alternate models, including a hierarchy of 
model fidelity. Although initially implemented with low- 
fidelity models, typical of the conceptual design 
environment, ultimately the architecture must allow 
multidisciplinary design, analysis, and optimization. 

Definition and development of the NDARC requirements 
benefited substantially from the experiences and computer 
codes of the preliminary design team of the U.S. Army 
Aeroflightdynamics Directorate (AFDD) at Ames 
Research Center. This background is described in 
reference 6. In the early 1990s, the RC code (for 
RotorCraft) emerged from the AFDD efforts, with RC97 a 
major version that unified tiltrotor and helicopter analyses. 
NDARC is entirely new software, built on a new 
architecture for the design and analysis of rotorcraft. From 
the RC theoretical basis, the equations of the parametric 
weight equations and the Referred Parameter Turboshaft 
Engine Model were used with only minor changes. Use 
was also made of some RC component aerodynamic 
models and rotor performance models. The current users 
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of RC, informed by past and recent applications, 
contributed significantly to the requirements definition. 

NDARC TASKS 

The NDARC code performs design and analysis tasks. 
The design task sizes the rotorcraft to satisfy specified 
design conditions and missions. The analysis tasks can 
include off-design mission performance analysis, flight 
performance calculation for point operating conditions, 
and generation of subsystem or component performance 
maps. Figure 1 illustrates the tasks. The principal tasks 
(sizing, mission analysis, flight performance analysis) are 
shown in the figure as boxes with heavy borders. Black 
arrows show control of subordinate tasks. 


fixed model or 
previous job or 

DESIGN previous case ANALYZE 



Figure 1. Outline of NDARC tasks. 

The aircraft description (figure 1) consists of all the 
information, input and derived, that defines the aircraft. 
This information can be the result of the sizing task; can 
come entirely from input, for a fixed model; or can come 
from the sizing task in a previous case or previous job. 
The aircraft description information is available to all 
tasks and all solutions (indicated by green arrows). 

The sizing task determines the dimensions, power, and 
weight of a rotorcraft that can perform a specified set of 
design conditions and missions. The aircraft size is 
characterized by parameters such as design gross weight, 
weight empty, rotor radius, and engine power available. 
The relationships between dimensions, power, and weight 


generally require an iterative solution. From the design 
flight conditions and missions, the task can determine the 
total engine power or the rotor radius (or both power and 
radius can be fixed), as well as the design gross weight, 
maximum takeoff weight, drive system torque limit, and 
fuel tank capacity. For each propulsion group, the engine 
power or the rotor radius can be sized. 

Missions are defined for the sizing task and for the 
mission performance analysis. A mission consists of a 
specified number of mission segments, for which time, 
distance, and fuel bum are evaluated. For the sizing task, 
certain missions are designated for design gross weight 
calculations, for transmission sizing, and for fuel tank 
sizing. The engines are sized to meet the power 
requirement of each segment. The mission parameters 
include mission takeoff gross weight and useful load. For 
specified takeoff fuel weight with adjustable segments, the 
mission time or distance is adjusted so the fuel required 
for the mission (burned plus reserve) equals the takeoff 
fuel weight. The mission iteration is on time or distance (if 
adjustable), or on fuel weight. 

Flight conditions are specified for the sizing task and for 
the flight performance analysis. For the sizing task, certain 
flight conditions are designated for design gross weight 
calculations, for transmission sizing, for maximum takeoff 
weight calculations, and for antitorque or auxiliary-thrust 
rotor sizing. The engines are sized to meet the power 
requirement of each condition. The flight condition 
parameters include gross weight and useful load. 

For flight conditions and mission takeoff, the gross weight 
can be maximized such that the power required equals the 
power available. 

A flight state is defined for each mission segment and each 
flight condition. The aircraft performance can be analyzed 
for the specified state, or a maximum- effort performance 
can be identified. The maximum effort is specified by 
identifying a target (such as best endurance, best range, or 
power required equal power available) to be achieved by 
adjusting a variable (such as speed, rate of climb, or 
altitude). 

The aircraft must be trimmed by solving for the controls 
and motion that produce equilibrium in the specified flight 
state. Different trim solution definitions are required for 
various flight states, hence for each mission segment and 
flight condition. 

Evaluating the rotor hub forces and blade pitch angles may 
require solution of the blade flap equations of motion. 
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The following sections describe the NDARC tasks in more 
detail. 

AIRCRAFT DESCRIPTION 

Decomposition of the aircraft system into fundamental 
components is critical to achieving the ability to rapidly 
model a wide array of rotorcraft concepts. Thus the 
aircraft consists of a set of components, including 
fuselage, rotors, wings, tails, and propulsion. For each 
component, attributes such as performance, drag, and 
weight can be calculated. The aircraft attributes are 
obtained from the sum of the component attributes. 
Description and analysis of conventional rotorcraft 
configurations is facilitated, while retaining the capability 
to model novel and advanced concepts. Specific rotorcraft 
configurations considered are single main-rotor and tail- 
rotor helicopter; tandem helicopter; coaxial helicopter; and 
tiltrotors. Novel and advanced concepts typically are 
modeled by starting with one of these conventional 
configurations. For example, compound rotorcraft can be 
constructed by adding wings and propellers. 

The following components form the aircraft. 

a) Systems: The systems component contains weight 
information (fixed useful load, vibration, contingency, 
systems and equipment) for the aircraft. 

b) Fuselage: There is one fuselage for the aircraft. 

c) Landing Gear: There is one landing gear for the aircraft. 

d) Rotors: The aircraft can have one or more rotors, or no 
rotors. In addition to main rotors, the component can 
model tail rotors, propellers, proprotors, and ducted fans. 
A rotor is designated a main rotor, tail rotor, or propeller; 
and can be tilting, ducted, and/or antitorque. The rotor 
power required is evaluated using the energy method, as a 
sum of induced, profile, and parasite power. The power 
components are calculated in terms of an induced power 
factor and a mean drag coefficient. The power models 
account for the influence of speed, thrust, compressibility, 
stall, lift offset, and the induced interference between twin 
rotors. The calibration of these induced and profile power 
models reflects the level of technology being considered. 
Blade element theory is used to calculate rotor hub forces 
and moments and to solve for blade pitch control or 
flapping. 

e) Forces: The force component is a simple model for a 
lift, propulsion, or control subsystem, including a weight 
and fuel flow description. 


f) Wings: The aircraft can have one or more wings, or no 
wings. 

g) Tails: The aircraft can have one or more horizontal or 
vertical tail surfaces, or no tails. 

h) Fuel Tank: There is one fuel tank component for the 
aircraft. There are one or more sizes of auxiliary fuel 
tanks. 

i) Propulsion Groups: There are one or more propulsion 
groups. Each propulsion group is a set of rotors and engine 
groups, connected by a drive system. The components 
define the power required, and the engine groups define 
the power available. There are one or more drive states, 
with a set of gear ratios for each state. The power required 
equals the sum of component power, transmission losses, 
and accessory losses. There are drive system torque limits, 
and rotor and engine shaft ratings. 

j) Engine Groups: Each engine group has one or more 
engines of the same type. For each engine type an engine 
model is defined. 

The engine model describes a particular engine, used in 
one or more engine groups. The engine performance 
information includes mass flow, fuel flow, jet thrust, and 
momentum drag at the required power. A Referred 
Parameter Turboshaft Engine Model (RPTEM) enables the 
aircraft performance analysis to cover the entire spectrum 
of operation. This model uses curve fits of referred 
performance from an engine deck, including the effect of 
turbine speed. The effects of size (a scaling model, based 
on mass flow) and technology (specific power and specific 
fuel consumption) are included in the engine model. 

The ability to define the aircraft control structure through 
input is a key feature for configuration generality. Aircraft 
controls are connected to component controls. Aircraft 
controls consist of pilot’s controls (for trim), configuration 
variables (e.g. tilt of nacelle/pylon, engine, rotor shaft), 
and direct connections to component controls. There can 
be one or more control states, with different connections to 
components (for example, to model the controls of a 
tiltrotor in helicopter mode and airplane mode flight). 
There are default control connections for each 
configuration. 

Weights are calculated or input for all components and 
subsystems. Calculated weights are obtained from 
parametric equations based on weights of existing turbine 
powered helicopters and tiltrotors (and some fixed wing 
aircraft component weights). Multiplicative technology 
factors can be used for the weights of all elements of the 
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aircraft, to either match measured weights or account for 
advanced technology. 

SIZING TASK 

The sizing task determines the dimensions, power, and 
weight of a rotorcraft that can perform a specified set of 
design conditions and missions. The aircraft size is 
characterized by parameters such as design gross weight 
W D or weight empty W E , rotor radius R , and engine 
power available P eng . The relationships between 
dimensions, power, and weight generally require an 
iterative solution. From the design flight conditions and 
missions, the task can determine the total engine power or 
the rotor radius (or both power and radius can be fixed), as 
well as the design gross weight, maximum takeoff weight, 
drive system torque limit, and fuel tank capacity. For each 
propulsion group, the engine power or the rotor radius can 
be sized. Alternatively, P cng and R can be input rather 
than sized. Several aircraft parameters can be determined 
by a subset of the design conditions and missions: 

a) Design gross weight W D . 

b) Maximum takeoff gross weight W^q . 

c) Drive system torque limit P DS limit . 

d) Fuel tank capacity W rucl _ cap . 

e) Antitorque or auxiliary-thrust rotor design thrust 
T 

± design * 

Alternatively, these parameters can be fixed at input 
values. The design gross weight W D can be fixed. The 
weight empty can be fixed (achieved by adjusting the 
contingency weight). 

For each flight condition and for each mission, the gross 
weight and useful load are specified. The gross weight can 
be input, or maximized, or fallout. For flight conditions, 
the payload or fuel weight can be specified, and the other 
calculated; or both payload and fuel weight specified, with 
gross weight fallout. For missions, the payload or fuel 
weight can be specified, the other fallout, and then time or 
distance of mission segments adjusted; or fuel weight 
calculated from mission, and payload fallout; or both 
payload and fuel weight specified (or payload specified 
and fuel weight calculated from mission), with gross 
weight fallout. 

Component Sizing 

Numerous choices are implemented for identification of 
independent (fixed) and dependent (fallout) design 
parameters, with parameter variation facilitated by 
automating dependencies. 


Engine Power. The engine size is described by the power 
P eng , which is the sea- level static power available per 
engine at a specified takeoff rating. 

Main Rotor. The main rotor size is defined by the radius 
R or disk loading WIA, thrust- weighted solidity a, hover 
tip speed V^ p , and blade loading C w /o = WlpAV* p o. 
With more than one main rotor, the disk loading and blade 
loading are obtained from an input fraction of design gross 
weight, W =f w W D . 

If the rotor radius is fixed for the sizing task, three of R 
(or WIA), C w la, V tip , and a are input; and the other 
parameters are derived. Optionally the radius can be 
calculated from a specified ratio to the radius of another 
rotor. 

If the sizing task determines the rotor radius (R and 
WIA), then two of C w !a, V tip , and a are input; and the 
other parameter is derived. The radius can be sized for just 
a subset of the rotors, with fixed radius for the others. 

Antitorque or Auxiliary-thrust Rotor. For antitorque and 
auxiliary-thrust rotors, three of R (or WIA), C w la, V tip , 
and a are input; and the other parameters are derived. 
Optionally the radius can be calculated from a specified 
ratio to the radius of another rotor. The disk loading and 
blade loading are based on sign , where 7^^ is the 
maximum thrust from designated design conditions. 
Optionally the tail rotor radius can be scaled with the main 
rotor radius and the disk loading. 

Wing. The wing size is defined by the wing area S or 
wing loading WIS; span (perhaps calculated from other 
geometry), chord, and aspect ratio. With more than one 
wing, the wing loading is obtained from an input fraction 
of design gross weight, W = f w W D . 

Two of the following parameters are input: area (or wing 
loading), span, chord, and aspect ratio; the other 
parameters are derived. The span can be calculated from 
the rotor radius, fuselage width, and clearance (typically 
used for tiltrotors). The span can be calculated from a 
specified ratio to the span of another wing. 

Tail. The tail size is determined by the area or tail volume, 
span, chord, and aspect ratio. 

Fuel Tank. The fuel tank capacity W fud-cap (maximum 
usable fuel weight) is determined from designated sizing 
missions. 

Weights. The structural design gross weight W SD and 
maximum takeoff weight W^q can be input, or specified 
as d + jW , for input increment d and fraction / . This 
convention allows the weights to be input directly ( / = 0), 
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or scaled with W D . For W SD , W is the design gross 
weight W D , or W D adjusted for a specified fuel state 
(input fraction of fuel capacity). For W^q, W is the 
design gross weight W D , or W D adjusted for maximum 
fuel capacity. Alternatively, can be calculated as 

the maximum gross weight possible at a designated sizing 
flight condition. 

Drive System Rating. The drive system rating is defined as 
a power limit, Poshmir The rating is properly a torque 
limit, Qosiumt = limit but is expressed as a power 
limit for clarity. The drive system rating can be input; 
calculated from the engine takeoff power rating; 
calculated from the power available at the transmission 
sizing flight conditions; or calculated from the power 
required at the transmission sizing flight conditions. The 
drive system rating is a limit on the entire propulsion 
system. To account for differences in the distribution of 
power through the drive system, limits are also used for 
the torque of each rotor shaft and of each engine group. 

MISSION DEFINITION 

Missions are defined for the sizing task and for the 
mission performance analysis. A mission consists of a 
specified number of mission segments. The takeoff gross 
weight is evaluated at the start of the mission, perhaps 
maximized for zero power margin at a specified mission 
segment. Then the aircraft is flown for all segments. For 
calculated mission fuel weight, the fuel weight at takeoff 
is adjusted to equal the fuel required for the mission 
(burned plus reserve). For specified takeoff fuel weight 
with adjustable segments, the mission time or distance is 
adjusted so the fuel required for the mission (burned plus 
reserve) equals the takeoff fuel weight. The mission 
iteration is thus on mission fuel weight. A successive 
substitution method is used if an iteration is required. 

Each mission consists of a specified number of mission 
segments. The following segment types can be specified. 

a) T axi or warm-up (fuel burned but no distance added to 
range). 

b) Distance: fly segment for specified distance (calculate 
time). 

c) Time: fly segment for specified time (calculate 
distance). 

d) Hold: fly segment for specified time (loiter, so fuel 
burned but no distance added to range). 

e) Climb: climb or descend from present altitude to next 
segment altitude (calculate time and distance). 

f) Spiral: climb or descend from present altitude to next 


segment altitude, fuel burned but no distance added to 
range. 

The number of auxiliary fuel tanks can change with each 
mission segment. The aircraft can refuel (either on the 
ground or in the air) at the start of a mission segment, by 
either filling all tanks to capacity or adding a specified fuel 
weight. Fuel can be dropped at the start of a mission 
segment. For calculation of the time or distance in a 
mission segment, a headwind or tailwind can be specified. 

Mission fuel reserves can be specified in several ways for 
each mission. Fuel reserves can be defined in terms of 
specific mission segments. Fuel reserves can be an input 
fraction of the fuel burned by all (except reserve) mission 
segments, or an input fraction of the fuel capacity. 

The takeoff distance can be calculated, either as ground 
run plus climb to clear an obstacle, or accelerate- stop 
distance in case of engine failure. Landing and VTOL 
takeoff calculations are not implemented, as these are best 
solved as an optimal control problem. 

FLIGHT STATE 

A flight state is defined for each flight condition and for 
each mission segment. The flight state definition consists 
of the speed, aircraft motion, altitude, atmosphere, height 
above ground level and landing gear state, aircraft control 
state, aircraft control values, and center-of-gravity 
position. Parameters defined for each propulsion group 
include drive system state, and rotor tip speed for primary 
rotor. Specified for each engine group are the number of 
inoperative engines, the infrared suppressor state, the 
engine rating, and the fraction of rated engine power 
available. Aircraft and rotor performance parameters for 
each flight state include payload drag, rotor performance, 
and aircraft trim state and trim targets. 

The aircraft performance can be analyzed for the specified 
state, or a maximum-effort performance can be identified. 
The maximum effort is specified by identifying a target 
(such as best endurance, best range, or power required 
equal power available) to be achieved by adjusting a 
variable (such as speed, rate of climb, or altitude). Two 
maximum effort quantity/variable pairs can be specified, 
solved in nested iterations. 

Environment and Atmosphere 

The aerodynamic environment is defined by the speed of 
sound c s , and density p, and kinematic viscosity 
v = p f p of the air (or other fluid). These quantities can be 
obtained from the standard day (International Standard 
Atmosphere), or input directly. The International Standard 
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Atmosphere (ISA) is a model for the variation with 
altitude of pressure, temperature, density, and viscosity, 
published as International Standard ISO 2533 by the 
International Organization for Standardization (ISO) (ref. 
8). The ISA consists of a series of altitude ranges with 
constant lapse rate (linear temperature change with 
altitude). 


Sizing Task 

Size Iteration 


method: successive 
substitution 


Flight Conditions 


Missions 


Mission Analysis 


Missions 


Flight Performance Analysis 


Flight Conditions 


Flight State 
Maximum Effort 


method: golden section search for maximum 
endurance, range, or climb; 
otherwise secant or false position 


Trim 

method: Newton-Raphson 

Component Performance Evaluation 

Blade Flapping 

method: Newton-Raphson 


Flight Condition 
Maximum GW 

method: secant or false position 


Flight State 


Mission 

Mission Iteration 

fuel weight, adjust time/distance 

method: successive substitution 
Segments 
Maximum GW 

method: secant or false 
position 

Flight State 


Figure 2. Design and analysis tasks, with nested loops and 
solution methods. 


SOLUTION PROCEDURE 

The NDARC code performs design and analysis tasks. 
Figure 1 illustrates the tasks. The principal tasks (sizing, 
mission analysis, flight performance analysis) are shown 
in the figure as boxes with dark borders. Black arrows 
show control of subordinate tasks. The aircraft description 
consists of all the information, input and derived, that 
defines the aircraft. This information can be the result of 
the sizing task; can come entirely from input, for a fixed 
model; or can come from the sizing task in a previous case 
or previous job. The aircraft description information is 


available to all tasks and all solutions (indicated by green 
arrows in figure 1). 

The nested iteration loops involved in the solution process 
are indicated by the subtitles in the boxes of figure 1 , and 
illustrated in more detail in figure 2. The flight state 
solution involves up to three loops. The innermost loop is 
the solution of the blade flap equations of motion, needed 
for an accurate evaluation of the rotor hub forces. The next 
loop is the trim solution, which is required for most flight 
states. The flight state optionally has one or two maximum 
effort iterations. The flight state solution is executed for 
each flight condition and for each mission segment. A 
flight condition solution or any mission segment solution 
can optionally maximize the aircraft gross weight. The 
mission usually requires an iterative solution, for fuel 
weight or for adjustable segment time or distance. Thus 
each flight condition solution involves up to four nested 
iterations: maximum gross weight (outer), maximum 
effort, trim, and blade motion (inner). Each mission 
solution involves up to five nested iterations: mission 
(outer), then for each segment maximum gross weight, 
maximum effort, trim, and blade motion (inner). Finally, 
the design task introduces a sizing iteration, which is the 
outermost loop of the process. The solution method for 
each iteration is indicated in figure 2. Details of the 
solution methods and their implementation are given in 
reference 6. 


COST 

Costs are estimated using statistical models based on 
historical aircraft price and maintenance cost data, with 
appropriate factors to account for technology impact and 
inflation. The aircraft flyaway cost consists of airframe, 
mission equipment package, and flight control electronics 
costs. The direct operating cost plus interest (DOC+I, in 
cents per available seat mile) is the sum of maintenance 
cost, flight crew salary and expenses, fuel and oil cost, 
depreciation, insurance cost, and finance cost. Inflation 
factors can be input, or internal factors used (either DoD 
deflators or CPI factors). 

The CTM rotorcraft cost model (refs. 9 to 11) gives an 
estimate of aircraft flyaway cost and direct operating cost 
plus interest. The statistical equation for airframe cost 
predicts the price of 123 out of 128 rotorcraft within 20% 
(figure 3). The fuel bum, block time, and block range are 
obtained from a designated mission. 
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Figure 3. Statistical estimation of rotorcraft flyaway cost. 


AIRCRAFT MODEL 

The aircraft disk loading is the ratio of the design gross 
weight and a reference rotor area: DL = W d / A ref . The 
reference area is a sum of specified fractions of the rotor 
areas, A rcf = 2 / a A. Typically A ref is the projected area 
of the lifting rotors, where the lifting rotors are all those 
not designated antitorque or auxiliary-thrust. The aircraft 
wing loading is the ratio of the design gross weight and a 
reference wing area: WL = W D fS K f . The reference area is 
a sum of the wing areas, *$ref “ 2 S . 

Aircraft Controls 

A set of aircraft controls c AC are defined, and connected 
to the component controls. The connection to the 
component control c is typically of the form 
c = STc ac + c 0 , where T is an input matrix and c 0 the 
component control for zero aircraft control. The factor S 
is available for internal scaling of the matrix. Figure 4 
illustrates the control relationships. 

The connection (matrix T) is defined for several control 
system states, allowing change of control configuration 
with flight state. An example is control of a tiltrotor in 
helicopter mode and airplane mode flight. The control 
state and initial control values are specified for each flight 
state. Default control matrices are defined based on the 
configuration: helicopter, tandem, coaxial, or tiltrotor. 


aircraft component controls 

controls c=Tc ac +c 0 

c AC c 0 = zero, constant, f (V) 



Figure 4. Aircraft and component controls 

Typical (default) aircraft controls are the pilot's controls: 
collective stick, lateral and longitudinal cyclic sticks, 
pedal, and tilt. Units and sign convention of the pilot's 
controls are contained in the matrix T . 

These aircraft controls are available for trim of the aircraft. 
Any aircraft controls not selected for trim will remain 
fixed at the values specified for the flight state. Thus by 
defining additional aircraft controls, component controls 
can be specified as required for a flight state. 

The tilt control variable a tilt is intended for nacelle tilt 
angle or conversion control, particularly for tiltrotor s. 
Typically this variable will be connected to nacelle/pylon, 
engine, rotor shaft, and even wings. The convention is 
a m = 0 for cruise, and a m = 90 deg for helicopter mode. 

Trim 

The aircraft trim operation solves for the controls and 
motion that produce equilibrium in the specified flight 
state. In steady flight (including hover, level flight, climb 
and descent, and turns), equilibrium implies zero net force 
and moment on the aircraft. There can be additional 
quantities that at equilibrium must equal target values. In 
practice, the trim solution can deal with a subset of these 
quantities. Usually it is at least necessary to achieve 
equilibrium in the aircraft lift and drag forces, as well as in 
yaw moment for torque balance. The basic purpose of the 
trim solution is to determine the component states, 
including aircraft drag and rotor thrust, sufficient to 
evaluate the aircraft performance. 

Different trim solution definitions are required for various 
flight states. Therefore one or more trim states are defined 
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for the analysis, and the appropriate trim state selected for 
each flight state of a performance condition or mission 
segment. For each trim state, the trim quantities, trim 
variables, and targets are specified. The available trim 
quantities include aircraft total force and moment, load 
factor, power or power margin, rotor force or torque, rotor 
flapping, wing force or lift margin, and tail force. The 
available trim variables include aircraft controls, 
orientation, speed, pullup rate, and turn rate. 

Geometry 

Layout of the geometry is typically in terms of station line 
(SL, positive aft), buttline (BL, positive right), and 
waterline (WL, positive up), measured relative to some 
arbitrary origin. The x, y, and z axes are parallel to the 
SL, BL, and WL directions. One or more locations are 
defined for each component of the aircraft. Each 
component will at least have a location that is the point 
where component forces and moments act on the aircraft. 
Each location is input in fixed or scaled form. The fixed 
form input is SL/BL/WL (dimensional). The scaled form 
input is jc /L (positive aft), y/L (positive right), and z/L 
(positive up), from a reference point. The reference length 
L is the rotor radius or wing span of a designated 
component, or the fuselage length. Fixed input can be used 
for the entire aircraft, or just for certain components. 

Aircraft Motion 

The aircraft velocity and orientation are defined by the 
flight speed, turn rate, orientation of the body frame 
relative to inertial axes (Euler angles), and orientation of 
the velocity frame relative to inertial axes (flight path 
angles). 

Aircraft conventions are followed for the direction and 
orientation of axes: the z-axis is down, the jc-axis 
forward, and the y-axis to the right; and a yaw-pitch-roll 
sequence is used for the Euler angles. 

The orientation of the body frame F relative to inertial 
axes is defined by yaw ip F , pitch d F , and roll <p F Euler 
angles. The flight path is specified by the velocity V, in 
the positive jc-axis direction of the velocity axes. The 
orientation of the velocity axes relative to inertial axes is 
defined by yaw ip v (sideslip) and pitch d v (climb) angles. 
In straight flight, all these angles are constant. In turning 
flight at a constant yaw rate, the yaw angle is ip F = ip F t . 
The aircraft angular velocity a) F c is obtained from the 
Euler angle rates. For steady state flight, 6 F =(j> F = 0; ip F 
is nonzero in a turn. Accelerated flight is also considered, 
in terms of linear acceleration and pitch rate. 


Loads and Performance 

For each component, the power required and the net forces 
and moments acting on the aircraft can be calculated. The 
aerodynamic forces F and moments M are typically 
calculated in wind axes and then resolved into body axes, 
relative to the origin of the body axes (the aircraft center- 
of-gravity). The power and loads of all components are 
summed to obtain the aircraft power and loads. Typically 
the trim solution drives the net forces and moments on the 
aircraft to zero. 

The aircraft equations of motion, in body axes F with 
origin at the aircraft center-of-gravity, are the equations of 
force and moment equilibrium. The forces and moments 
are the sum of loads from all components of the aircraft: 
fuselage, rotor, force, wing, tail, engine, and tank. A 
particular component can have more than one source of 
loads; for example, the rotor component produces hub 
forces and moments, but also includes hub and pylon drag. 

The component power required P comp is evaluated for all 
components (rotors) of the propulsion group. The total 
power required for the propulsion group P req pQ is 
obtained by adding the transmission losses and accessory 
power. The power required for the propulsion group must 
be distributed to the engine groups. From the power 
required the fuel flow is calculated. The fuel flow of the 
propulsion group is obtained from the sum over the engine 
groups. The total fuel flow is the sum from all components 
of the aircraft: w = Zw reqEG + 2w force 

Aerodynamics 

Each component has a position z F in aircraft axes F, 
relative to the reference point, and orientation of 
component axes B relative to aircraft axes given by a 
rotation matrix. Acting at the component are interference 
velocities (velocity of air, in F axes) from all other 
components. Then the total component velocity relative to 
the air is the sum of velocity, angular velocity, and 
interference terms. The aerodynamic environment is 
defined in the component axes: velocity magnitude 
v =1 v B I, dynamic pressure q = %pv 2 , angle-of-attack a, 
and sideslip angle /L From these the component 
aerodynamic model calculates the loads in wind axes: 
drag, side force, and lift ( D, Y, and L ) and the roll, pitch, 
and yaw moments; and from wind axis loads the 
aerodynamic loads in aircraft axes acting at the center-of- 
gravity are calculated. 

The component aerodynamic models are based on 
equations intended to cover the principal operating 
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conditions, including small angle-of-attack and sideslip 
with a representation of stall and maximum lift, vertical 
flight, sideward flight, and rearward flight. 

Drag 

Each component can contribute drag to the aircraft. A 
fixed drag can be specified, as a drag area Dfq\ or the 
drag can be scaled, specified as a drag coefficient Q? 
based on an appropriate area S. There may also be other 
ways to define a scaled drag value. For fixed drag, the 
coefficient is C D ={Dlq)IS (the aerodynamic model is 
formulated in terms of drag coefficient). For scaled drag, 
the drag area is D!q = SC D . For all components, the drag 
(Plq ) or C Dcomp is defined for forward flight or 
cruise; typically this is the minimum drag value. For some 
components, the vertical drag or sideward drag is defined. 
For some components, the aerodynamic model includes 
drag due to lift, angle-of-attack, or stall. 

Table 1 summarizes the component contributions to drag, 
and the corresponding reference areas. If no reference area 
is indicated, then the input is only drag area Dlq. An 
appropriate drag reference area is defined for each 
component; the reference area is either input or calculated. 

Table 1. Component contributions to drag. 


component 

drag contribution 

reference area 

fuselage 

fuselage 

fuselage wetted area 


fuselage vertical 

fuselage projected area 


fixtures and 

fuselage wetted area 


fittings 

rotor-body interf 

fuselage wetted area 


contingency 

— 


payload 

increment 

— 

gear 

landing gear 

— 

rotor 

hub, hub vertical 

rotor disk area 


pylon, pylon vert 

pylon wetted area 


spinner 

spinner wetted area 

wing 

wing, wing vert 

wing planform area 


wing-body interf 

wing planform area 

tail 

tail, tail vertical 

tail planform area 

engine 

nacelle, nac. vert 

nacelle wetted area 


momentum drag 

— 

fuel tank 

auxiliary tank 

— 


Optionally the aircraft drag can be fixed. The quantity 
specified is the sum (over all components) of the drag area 


Dlq (minimum drag, excluding drag due to lift and 
angle-of-attack), without accounting for interference 
effects on dynamic pressure. The input parameter can be 
D!q\ or the drag can be scaled, specified as a drag 
coefficient based on the rotor disk area, so D!q = A ref C D 
( A ref is the reference rotor disk area); or the drag can be 
estimated from the aircraft maximum take-off weight, 
Dlq = kfW^Q /1000) 2/3 . Based on historical data, the 
drag coefficient C D = 0.02 for old helicopters and 
C D = 0.008 for current low drag helicopters. Based on 
historical data, k = 9 for old helicopters, k = 2.5 for 
current low drag helicopters, £ = 1.6 for current tiltrotors, 
and k = 1.4 for turboprop aircraft (k in English units). If 
the aircraft drag is input, then the fuselage contingency 
drag is adjusted so the total aircraft Dlq equals the input 
value. 

The nominal drag areas of the components and the aircraft 
are part of the aircraft description and are used when the 
aircraft drag is fixed. While vertical drag parameters are 
part of the aerodynamic model for the hub, pylon, and 
nacelle, aerodynamic interference at the rotor and at the 
engine group is not considered, so these terms do not 
contribute to download. In the context of download, only 
the fuselage, wing, tail, and contingency contribute to the 
nominal vertical drag. 

FUSELAGE 

There is one fuselage component for the aircraft. There are 
several options for calculating the fuselage length, wetted 
area (reference area for drag coefficients), and projected 
area (reference area for vertical drag). 

The aerodynamic velocity of the fuselage relative to the 
air, including interference, is calculated in component 
axes. The drag area or drag coefficient is defined for 
forward flight, vertical flight, and sideward flight. For 
small angle-of-attack, the drag increases proportional to 
(1 + K\a\ x ), using input factor K and exponent X. In 
addition, the forward flight drag area or drag coefficient is 
defined for fixtures and fittings, and for rotor-body 
interference. The drag force is 

^ — *7*^wet + ^Dfit + ^^Drb) + #)pay + iP ! 4) cont) 

including the drag area of the payload (specified for the 
flight state) and the contingency drag area. 

The fuselage lift, pitch moment, side force, and yaw 
moment are defined in fixed form (e.g. Llq) or scaled 
form (e.g. C L ). Maximum lift coefficient and maximum 
side force coefficient are defined. 
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ROTOR 

The aircraft can have one or more rotors, or no rotors. In 
addition to main rotors, the rotor component can model 
tail rotors, propellers, proprotors, ducted fans, thrust 
vectoring rotors, and auxiliary-thrust rotors. The principal 
configuration designation (main rotor, tail rotor, or 
propeller) determines where the weights are put in the 
weight statement (rotor group, empennage group, or 
propulsion group), and each configuration can possibly 
have a separate performance or weight model. Antitorque 
rotors and auxiliary-thrust rotors can be identified, for 
special sizing options. Other configuration features are 
variable diameter and ducted fan. 

Multi-rotor systems (such as coaxial or tandem 
configuration) are modeled as a set of separate rotors, in 
order to accommodate the description of the position, 
orientation, controls, and loads. The performance 
calculation for twin rotor systems can include the mutual 
influence of the induced velocity on the power. 

Drive System 

The drive system defines gear ratios for all the 
components that it connects (rotors and engine groups). 
There is one primary rotor per propulsion group (for 
which the reference tip speed is specified); other 
components are dependent (for which a gear ratio is 
specified). For the primary rotor, a reference tip speed 
V ti p-ref is defined for each drive system state. 

Geometry 

The rotor rotation direction is described by the parameter 
r: r = l for counter-clockwise rotation, r = -l for 

clockwise rotation (as viewed from the positive thrust side 
of the rotor). The rotor solidity and blade mean chord are 
related by a = Nc I JtR . Generally thrust- weighted values 
are used, but geometric values are also required by the 
analysis. A general blade chord distribution is specified as 
c(r) = c T&f c(r), where c ref is the thrust- weighted chord. 
Linear taper is specified in terms of a taper ratio 
? = c tip /c root , or in terms of the ratio of thrust- weighted 
and geometric chords, / = a t fa g = c 75R fc 50 R . 

The rotor hub is at position z bub . A component of the 
position can be calculated, superseding the location input. 
The calculated geometry depends on the configuration. 
For a coaxial rotor, the hub locations can be calculated 
from the input separation, and the input location midway 
between the hubs. For a tandem rotor, the hub locations 
can be calculated from the input overlap, and the input 
location midway between the hubs. For a tail rotor, the 


longitudinal position can be calculated from the main rotor 
radius, tail rotor radius, and tail-rotor/main-rotor 
clearance. For a tiltrotor, the lateral position can be 
calculated from the rotor radius (cruise value for variable 
diameter rotor), fuselage/rotor clearance, and fuselage 
width (with the pivot, pylon, and nacelle center-of-gravity 
lateral positions adjusted to keep the same relative position 
to the hub). Alternatively for a tiltrotor, the lateral position 
can be calculated from the wing span so the rotors are at 
the wing tips; or from a designated wing panel edge. 

Control and Loads 

The rotor controls consist of collective, lateral cyclic, 
longitudinal cyclic, and perhaps shaft incidence (tilt) and 
cant angles. Rotor cyclic control can be defined in terms 
of tip-path plane or no-feathering plane command. The 
collective control variable is the rotor thrust amplitude T 
or C T la (in shaft axes), from which the collective pitch 
angle can be calculated. This approach eliminates an 
iteration between thrust and inflow. 

The relationship between tip-path plane tilt and hub 
moment is M = (N !2)I b Q 2 (v 2 -1)/? = K huh fi, where N 
is the number of blades, Q the rotor speed, and v the 
dimensionless fundamental flap frequency. The flap 
moment of inertia I b is obtained from the input Lock 
number. 

Optionally the rotor can have a variable diameter. The 
rotor diameter is treated as a component control. As a 
control it can be connected to an aircraft control and thus 
set for each flight state. 

Tip-path plane command is characterized by direct control 
of the rotor thrust magnitude and the tip-path plane tilt. 
Cyclic tilt of the tip-path plane, hence tilt of the thrust 
vector, consists of longitudinal tilt (positive forward) 
and lateral tilt f$ s (positive toward retreating side). 
Alternatively, the cyclic control can be specified in terms 
of hub moment or lift offset, if the blade flap frequency is 
greater than 1/rev. This control mode requires calculation 
of rotor collective and cyclic pitch angles from the thrust 
magnitude and flapping. 

No-feathering plane command is characterized by control 
of rotor cyclic pitch angles, and direct control of the rotor 
thrust magnitude. Cyclic tilt of the no-feathering plane, 
usually producing tilt of the thrust vector, consists of 
longitudinal cyclic pitch angle 6 S (positive aft) and lateral 
cyclic pitch angle 6 C (positive toward retreating side). 
This control mode requires calculation of rotor collective 
pitch angle and tip-path plane tilt from the thrust 
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magnitude and cyclic control, including the influence of 
inflow. 

Rotor Axes and Shaft Tilt 

The rotor hub is at position z£ ub , where the rotor forces 
and moments act; the orientation of the rotor shaft axes 
relative to the aircraft axes is given by the rotation matrix 
C SF . The pivot is at position Zp ivot . The hub or shaft axes 
S have origin at the hub center; the z-axis is the shaft, 
positive in the positive thrust direction; and the jc-axis 
usually downstream. The rotor orientation is specified by 
selecting a nominal direction in body axes (positive or 
negative x, y, or z-axis) for the positive thrust direction; 
the other two axes are then the axes of shaft control. For a 
main rotor the nominal direction would be the negative z- 
axis; for a tail rotor the nominal direction would be the 
±y-axis (depending on the direction of rotation of the 
main rotor); for a propeller the nominal direction would be 
the positive x-axis. The hub and pivot axes have a fixed 
orientation relative to the body axes: hub incidence and 
cant, pivot dihedral, pitch, and sweep. The shaft angle 
control consists of incidence and cant about the pivot axes. 

For a tiltrotor aircraft, one of the aircraft controls is the 
nacelle angle, with the convention a m = 0 for cruise, and 
a m = 90 deg for helicopter mode. The rotor shaft 
incidence angle is then connected to a m by defining the 
control matrix appropriately. 

Hub Loads 

The rotor controls give the thrust magnitude and the tip- 
path plane tilt angles and , either directly or from the 
collective and cyclic pitch. The forces acting on the hub 
are the thrust T , drag H , and side force Y (positive in z, 
x, and y-axis directions respectively). The hub pitch and 
roll moment are proportional to the flap angles. The hub 
torque is obtained from the shaft power P shaft and rotor 
speed Q. The force and moment acting on the hub, in 
shaft axes, are then: 


(H\ 


° ' 

Y 

+ 

0 

Ji 


-/bTj 


'M x \ (K hub (rp s )> 

m s = M y = tf hub (-/y 

\~ r Q/ \ ^shaft 1^) 


drag of the fuselage or a lifting surface in the aircraft trim. 
For example, f B can model the tail rotor blockage caused 
by operation near the vertical tail. 

The rotor loads in aircraft axes acting at the center-of- 
gravity ( F f and M F ) are then calculated from the shaft 
axis loads ( F s and M s ) and Az = Zhub “ z cg • 

The wind axes lift L and drag X are calculated from the 
net rotor hub force F f and the rotor velocity v F . The 
velocity relative to the air gives the propulsive force 
direction e P =v F l\v F \ (no interference), and the 
velocity magnitude V =1 v F I. The lift and drag 
components of the force are XV = -(v F ) T F f and 
L 2 =1 F f I 2 -X 2 . 

Aerodynamics 

The rotor velocity relative to the air is 
v F =v ac +(0 ac *Az F in aircraft axes. The velocities in 
shaft axes are 

(-*) 

V s =C SF v F =QR rn y 

j 

where QR is the rotor tip speed. The advance ratio, inflow 
ratio, and shaft angle-of-attack are defined as 
^ + tiy) V2 , A = A f + j u z , and a = tan _1 (ju z /ju). The 

rotor thrust coefficient is defined as C T =T / pA(QR ) 2 . 
The dimensionless ideal induced velocity A f is calculated 
from ju, \i v and C T ; the dimensional velocity is 
v ( =QRh i . The ideal induced power is then P ideal =Tv ( . 
Note that for these inflow velocities, the subscript “ i ” 
denotes “ideal.” 

Ideal Inflow 

The ideal wake-induced velocity is obtained from the 
momentum theory result of Glauert: 



where A = A* + ju z , A^ =1 C T I /2 ( X h is always positive), 
and s = sign(C r ). This expression is generalized to 
Aj* = X h sF(ii! k hi sii z ! k h ). If jU is zero, the equation for 
A; can be solved analytically. Otherwise, for non-axial 
flow, an iterative Newton-Raphson solution for A f is 
needed. 


The force includes a term equal to the rotor thrust times an 
input blockage factor f B = AT IT . This term accounts for 
blockage or download, as an alternative to including the 
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An approximate equation is used in the turbulent-wake 
and vortex-ring states to eliminate the singularity of the 
momentum theory result at ideal autorotation. The 



momentum theory result is also extended to the case of a 
ducted fan. 

The wake-induced velocity is reduced when the rotor disk 
is in the proximity of the ground plane, 
(^j)ige = f g (^j)oge* The factor / is a function of the 
scaled rotor height above the ground, z g ID. The effects 
of ground plane tilt and rotor velocity are modeled. 
Several empirical ground effect models are implemented. 

As a simple approximation to nonuniform induced 
velocity distribution, a linear variation over the disk is 
used: AA = A x r cos xp + k y r sin rp . There are contributions 
to AA from forward flight and from hub moments, which 
influence the relationship between flapping and cyclic. 
The models implemented for the forward flight gradients 
are based on references 12 to 15. Differential momentum 
theory is used to calculate the gradients caused by hub 
moments. 


significant, as for a rotor with large flap stiffness. Figure 5 
shows the tip-path plane tilt and thrust vector tilt with 
cyclic pitch control (no-feathering plane tilt) as a function 
of flap stiffness (per-rev frequency) for several rotor thrust 
values. The difference between tip-path plane tilt and 
thrust vector tilt is caused by tilt of the thrust vector 
relative to the tip-path plane. 

The profile inplane forces can be obtained from simplified 
equations, or calculated by blade element theory. The 
simplified method uses: 

Cho\_(^ p 

Cyj-r*™ 

where the mean drag coefficient c dmem is from the profile 
power calculation. The function F„ accounts for the 
increase of the blade section velocity with rotor edgewise 
and axial speed. 


Rotor Forces 

Direct control of the rotor thrust magnitude is used, so the 
rotor collective pitch angle 0 O must be calculated from the 
thrust C T la . If the commanded variable were the 
collective pitch angle, then calculating the rotor thrust 
would be necessary, resulting in a more complicated 
solution procedure; in particular, an iteration between 
thrust and inflow would be needed. There may be flight 
states where the commanded thrust can not be produced 
by the rotor, even with stall neglected in the section 
aerodynamics. This situation will manifest as an inability 
to solve for the collective pitch given the thrust. In this 
circumstance the trim method should be changed so the 
required or specified thrust is an achievable value. 

The inplane hub forces are produced by tilt of the thrust 
vector with the tip-path plane, plus forces in the tip-path 
plane, and profile terms (produced by the blade drag 
coefficient). The orientation of the tip-path axes relative to 
the shaft axes is C PS , which is determined by and fi s . 
Then 



The inplane forces relative to the tip-path plane can be 
neglected or calculated by blade element theory. Note that 
with tip-path plane command and C Htp and C Ftpp 
neglected, solving for the rotor collective and cyclic pitch 
angles is not necessary. In general the inplane forces 
relative to the tip-path plane are not zero, and may be 



flap frequency v (per-rev) 



1.0 1.2 1.4 1.6 1.8 2.0 


flap frequency v (per-rev) 

Figure 5. Tip-path plane tilt and thrust vector tilt with 
cyclic pitch. 
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Blade Element Theory 

Blade element theory is the basis for the solution for the 
collective and cyclic pitch angles (or flap angles) and 
evaluation of the rotor inplane hub forces. The section 
aerodynamics are described by lift varying linearly with 
angle-of-attack, c £ = c ia a (no stall), and a constant mean 
drag coefficient c rfmean (from the profile power 
calculation). The blade section aerodynamic environment 
is described by the three components of velocity, from 
which the yaw and inflow angles are obtained, and then 
the angle-of-attack a . The blade pitch consists of 
collective, cyclic, twist, and pitch-flap coupling terms. The 
flap motion is rigid rotation about a hinge with no offset, 
and only coning and once-per-revolution terms are 
considered. The inflow includes gradients caused by 
edgewise flight and hub moments. The effect of the inflow 
produced by hub moments is to introduce a lift- deficiency 
function in the flap response. 

Integrating the section lift, drag, and radial forces over the 
blade radial coordinate and azimuth gives the total rotor 
thrust, drag, and side forces. Integrating terms produced 
by the section drag coefficient gives the profile inplane 
forces, C Ho and C Yo \ using blade element theory to 
evaluate these accounts for the planform and root cutout, 
while using F„ implies a rectangular blade and no root 
cutout (plus at most a 1% error approximating the exact 
integration). The remaining terms in the section forces 
produce the inplane loads relative to the tip-path plane, 

^tftpp an d Crtpp* 

Evaluating these inplane forces requires the collective and 
cyclic pitch angles, and the flapping motion. The thrust 
equation must be solved for the rotor collective pitch. The 
relationship between cyclic pitch and flapping is defined 
by the rotor flap dynamics. The flap motion is rigid 
rotation about a central hinge, with a flap frequency v > 1 
for articulated or hingeless rotors. The flapping equation 
of motion, including the effects of precone and the inertial 
loads of shaft angular motion, is harmonically analyzed. 
The result is equations for the mean (coning) and 1/rev 
(tip-path plane tilt) flap motion. The solution for the 
coning is largely decoupled by introducing the rotor thrust 
coefficient. 

The thrust and flapping equations of motion must be 
solved for the unknown angles. For tip-path plane 
command, the thrust and flapping are known, so the 
equations are solved for d 0 75 , 6 C , and 6 S . For no- 
feathering plane command, the thrust and cyclic pitch are 
known, so the equations are solved for d 0 75 , and /3 S . 


Power 

The rotor power consists of induced, profile, and parasite 
terms: P = P t + P Q + P p . The rotor parasite power 

(including climb/descent power for the aircraft) is 
obtained from the wind axis drag force: 

P p =-XV = (v f ) t F f . 

The induced power is calculated from the ideal power: 
p. = KP idc?d = K/£>7v ideal . The empirical factor k accounts 
for the effects of nonuniform inflow, non-ideal span 
loading, tip losses, swirl, blockage, and other phenomenon 
that increase the induced power losses (ff>l). For a 
ducted fan, f D =f w ( 2 is introduced, where f w is the 
ratio of the far wake induced velocity to the induced 
velocity at the disk. 

The profile power is calculated from a mean blade drag 
coefficient: P a = pA(QR) 3 C Po , C Po = (a/8)c dmean F P . 

The function F p (ja,ia z ) accounts for the increase of the 
blade section velocity with rotor edgewise and axial speed. 

Two performance methods are implemented: the energy 
method and the table method. The induced power factor 
and mean blade drag coefficient are obtained from 
equations with the energy method, or from tables with the 
table method. Optionally k and c rfmean can be specified 
for each flight state, superseding the performance method 
values. 

Energy Performance Method: Induced Power 
The induced power is calculated from the ideal power: 
P i= KP id e ai= K fD T Vi< ie ai- Reference values of k are 
specified for hover, axial cruise (propeller), and edgewise 
cruise (helicopter): /c hovcr , p , k^. Two models are 
implemented: constant model and standard model. The 
constant model uses k = ff hover if jU = \i z = 0; or k = K prop 
if I jU k 0.1 1 \x z I; or k = otherwise. 

The standard model calculates an axial flow factor 
from the values K hover , K cUmb , and K pmp . Let 
A = C T la- (C T Icr)^. For hover and low speed axial 
climb, including a variation with thrust, the inflow factor 
is 

K h = hover + ^hl^h + ^h2^h 

"I" 0^* climb _ * hover)! tan' 1 [((I p z \/X h )/M 

axial) ] 

where I \i z I / k h = is the midpoint of the transition 

between hover and climb. Figure 6 illustrates k in hover 
(with a minimum value). A polynomial describes the 
variation with axial velocity, scaled so k = K h at ju z = 0 
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and K=K p at \i z = M zprop > including a variation with 
thrust: 


k p ~ ^prop + kpi^p + k p2 A p 
K axial = K h + K\ P Z + (KlPz + KsPz' ° ) 


where £ a accomplishes the scaling. A polynomial 
describes the variation with edgewise advance ratio, scaled 
so k= at ju = 0 and k = f off K e dge at \x = ju^. Thus 
the induced power factor is 

K = ^ axial + k elV + ( Kl ^ ) 

where accomplishes the scaling. The function 

/off = 1 - k ol (1 - e“* o20jc ) accounts for the influence of lift 
offset, o x = rM x ITR = (K huh /TR)P s . Figure 7 illustrates 
k: in edgewise flight. Minimum and maximum values of 
the induced power factor are also specified. 



Cj'/o 


Figure 6. Induced power factor in hover (with minimum 
k = 1 . 12 ). 



Figure 7. Induced power factor in edgewise flight. 


Energy Performance Method: Profile Power 

The profile power is calculated from a mean blade drag 
coefficient: C Po = (o78)c rfmean F P . Since the blade mean 
lift coefficient is c t & 6C r to , the drag coefficient is 


estimated as a function of blade loading C T la (using 
thrust-weighted solidity). With separate estimates of the 
basic, stall, and compressibility drag, the mean drag 
coefficient is 


C d mean X^( c d basic + C d stall + C d comp ) 

where % is a technology factor. The factor 
S = (Re mi ! Re) 02 accounts for Reynolds number effects 
on the drag coefficient; Re is based on the thrust- 
weighted chord, 0.751^, and the kinematic velocity of 
the flight state; and Re ref corresponds to the input drag 
coefficient information. Array and equation models are 
implemented for the basic drag. In the array model the 
basic drag c d basic is input as a function of C T fa (linearly 
interpolated array). 


In the equation model the basic drag c d basic is a quadratic 
function of C T !o, plus an additional term allowing faster 
growth at high (sub- stall) angles of attack. Let 
A =1 C T la - (Cy la) D min I , where (C T la) Dmin 
corresponds to the minimum drag; and 
Asep Ct 1° I ~(Pt / a )sep • Values of the basic drag 
equation are specified for helicopter (hover and edgewise) 
and propeller (axial climb and cruise) operation: 


c dh - ^Ohel + ^lhelA + ^hel^ 2 + ^ S epA S ep P 

C dp = ^Oprop + ^1 prop A + ^2 prop A + ^sepAsep P 


The separation term is present only if A scp >0. The 
helicopter and propeller values are interpolated as a 
function of i u z : 

C d basic =C dh + ( C dp _C dft)“| tan ^ Pz ^ ^h) 

so I \x z I / k h = 1 is the midpoint of the transition. 

The stall drag increment represents the rise of profile 
power caused by the occurrence of significant stall on the 
rotor disk. Let A s =1 C T la\-(f s lf oS )(C T fo) s , where f s 
is an input factor. The function f off = 1 - d ol (1 - e ~ d ° 2 ° x ) 
accounts for the influence of lift offset, 
o x = rM x ITR = (K hub /TR)P S . Then 

c d stall = ^ + 

(zero if A^ =s 0). The blade loading at which the stall 
affects the entire rotor power, (i C T la ) s , is an input 
function of the velocity ratio V = (j u 2 + ju 2 ) 1/2 . 

Figure 8 shows typical stall functions (i C T la) s for two 
rotors with different airfoils. For reference, typical rotor 
steady and transient load limits are also shown. 

The compressibility drag increment depends on the 
advancing tip Mach number M at . Drag divergence and 
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similarity models are implemented. For the drag 
divergence model, let AM = M at - M M , where M dd is 
the drag divergence Mach number of the tip section. Then 
the compressibility increment in the mean drag coefficient 
is 

c d co mp = d ml AM + d m2 AM x - 

(ref. 16). M dd is a function of the advancing tip lift 
coefficient (available from blade element theory) and the 
tip airfoil thickness-to-chord ratio. 

Figure 9 illustrates the mean drag coefficient in hover, 
showing c ^ without and with the separation term, and the 
total for the high-stall and low-stall cases. Figure 10 
illustrates the mean drag coefficient in forward flight, 
showing the compressibility term c dcomp , and the growth 
in profile power with C T / a and jU as the stall drag 
increment increases. 


transient limit 

steady limit 

— high stall 



Figure 8. Profile power stall function. 
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Figure 9. Mean drag coefficient in hover. 
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Figure 10. Mean drag coefficient in forward flight; (a) 
high stall; (b) low stall. 


Twin Rotors 

For twin rotors, the induced power is determined by the 
induced velocity of the rotor system, not the individual 
rotors. The induced power is still obtained using 
p i = ^idcai = ^/o^ideai for each rotor, but the ideal 
induced velocity is calculated for an equivalent thrust ^Te 
based on the thrust and geometry of both rotors (see refs. 4 
and 6). The profile power calculation is not changed for 
twin rotors. 

Interference 

The rotor can produce aerodynamic interference velocities 
at the other components (fuselage, wings, tails). The 
induced velocity at the rotor disk is kv-, acting opposite 
the thrust (z-axis of tip-path plane axes). So 
v£ d = ~k P Kv i , and v P nd = C FP v£ d . The total velocity of 
the rotor disk relative to the air consists of the aircraft 
velocity and the induced velocity from this rotor: 
v total= vF - v L- The direction of the wake axis is thus 
e p = -C PF v P otsd / 1 v** I . The angle of the wake axis from 
the thrust axis is x = cos -1 I (k p ) T e p I. 
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The interference velocity v£ t at each component is 
proportional to the induced velocity V ind (and is in the 
same direction), with factors accounting for the stage of 
wake development and the position of the component 
relative to the rotor wake: 

v int = ^int/w/z/r/f V ind 

The factors f w f z account for axial development of the 
wake velocity, the factor f r accounts for immersion in the 
wake, and A" int is an input empirical factor. An additional 
factor f t for twin rotors is included. 

Optionally the development along the wake axis can be a 
step function ( f w f z = 0, 1, or f w for above the rotor, on 
the rotor disk, or below the rotor disk, respectively); or the 
wake develops with the nominal or input rate of change. 
Optionally the wake immersion can use the contracted 
radius or the uncontracted radius; be a step function 
( f r = 1 and 0 inside and outside the wake boundary); be 
always immersed ( f r = 1 always); or use an input 
transition distance. 

The interference factor K mi can be reduced from an input 
value at low speed to zero at high speed, with linear 
variation over a specified speed range. To account for the 
extent of the wing or tail area immersed in the rotor wake, 
the interference velocity can be calculated at several points 
along the span and averaged. 

Drag 

The rotor component includes drag forces acting on the 
hub and spinner and on the pylon. 

The hub drag can be fixed, specified as a drag area Dlq\ 
scaled, specified as a drag coefficient Cd based on the 
rotor disk area; or estimated based on the gross weight, 
using either a squared-cubed relationship or a square-root 
relationship. Based on historical data, the drag coefficient 
C D = 0.004 for typical hubs, C D = 0.0024 for current low 
drag hubs, and C D = 0.0015 for faired hubs. For the 
squared-cubed relationship : (D / q) hub = kfW^Q / 1 000) 2 1 3 
where W^q is the maximum take-off gross weight. 
Based on historical data, k = 1.4 for typical hubs, k = 0.8 
for current low drag hubs, and k = 0.5 for faired hubs ( k 
in English units). For the square-root relationship: 
( D/ q)h»b=k(W M ro /N r 0 i 0 T ) U2 where W ^ /N mtor is the 
maximum take-off gross weight per lifting rotor. Based on 
historical data (ref. 17), k = 0.074 for single rotor 
helicopters, £ = 0.049 for tandem rotor helicopters 
(probably a blade number effect), k = 0.038 for hingeless 
rotors, and k = 0.027 for faired hubs (k in English units). 


The pylon forward flight drag and vertical drag are 
specified as drag area or drag coefficient, based on the 
pylon wetted area. The spinner drag is specified as drag 
area or drag coefficient, based on the spinner wetted area. 

FORCE 

The force component is an object that can generate a force 
acting on the aircraft, possibly used for lift, or propulsion, 
or control. The amplitude of the force can be a fixed value, 
or connected to an aircraft control for trim. The direction 
of the force can be fixed, or connected to aircraft control. 

The force generation requires a fuel flow that is calculated 
from an input thrust specific fuel consumption: 
w = A(sfc) , where A is the force amplitude, and the units 
of sfc are pound/hour/pound or kilogram/hour/Newton. 

The force component weight is identified as either engine 
system or propeller/fan installation weight, both of the 
propulsion group. The force component weight is 
calculated from a specific weight and the design maximum 
force F max , plus a fixed increment: W = SF max + AW . 

WING 

The aircraft can have one or more wings, or no wings. The 
wing is described by planform area S, span b , mean chord 
c = S lb, and aspect ratio AR = b 2 IS. These parameters 
are for the entire wing. The geometry is specified in terms 
of two of the following parameters: S or wing loading 
W IS, b (perhaps calculated from other geometry), c, 
AR. With more than one wing, the wing loading is 
obtained from an input fraction of design gross weight, 
W = f w W D . Optionally the span can be calculated from a 
specified ratio to the span of another wing. 

For the tiltrotor configuration, the wing span can be 
calculated from the fuselage and rotor geometry (rotor 
radius, rotor-fuselage clearance, and fuselage width). The 
wing span can be calculated from the rotor hub position 
(regardless of how the rotor position is determined). 
Optionally the wing span can be calculated from an 
appropriate specification of all wing panel widths. 

Panels 

The wing planform is defined in terms of one or more 
wing panels (figure 11). Symmetry of the wing is 
assumed. Each panel has a straight aerodynamic center 
and linear taper. The aerodynamic center locus (in wing 
axes) is defined by sweep, dihedral, and offsets at the 
inboard edge relative to the aerodynamic center of the 
previous panel. The wing position is the mean 
aerodynamic center. 
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A panel is characterized by span b p (each side), mean 
chord c p , and area S p = 2 b p c p (both sides). The taper is 
defined by inboard and outboard chord ratios. The span for 
each panel (if there are more than two panels) can be a 
fixed input, a fixed ratio of the wing span, or free. The 
panel outboard edge (except for the wing tip) can be at a 
fixed input position, at a fixed station, calculated from the 
fuselage and rotor geometry, calculated from the hub 
position, or adjusted. The specification of panel spans and 
panel edges must be consistent, and sufficient to determine 
the wing geometry. To complete the definition of the 
geometry, one of the following quantities is specified for 
each panel: panel area S p ; ratio of panel area to wing 
area, S p /S', panel mean chord c p ; ratio of panel mean 
chord to wing mean chord, c p !c ; inboard and outboard 
chord ratios; or free. The total wing area equals the sum of 
all panel areas. 



edge 

nip 
^ ip 
c ip 


*ip 

Z IP 


edge 

VEp 

Vop 

^ Op 
C Op 


outboard panel edge, rj = y/(b/2) 

wing station 

chord ratio, A = c/c ret 

chord 

sweep (+ aft) 
dihedral (+ up) 

aero center offset (inboard, + aft) 
aero center offset (inboard, + up) 


span (each side) 
mean chord 
area = 2b p c p 


Aerodynamics 

The wing vertical drag can be fixed, specified as a drag 
area ( Dlq ) v ; or scaled, specified as a drag coefficient 
C DV based on the wing area; or calculated from an airfoil 
section drag coefficient (for -90 deg angle-of-attack) and 
the wing area immersed in the rotor wake (including 
changes in wing area due to flap and flaperon deflection). 

The wing lift, pitch moment, and roll moment are defined 
in scaled form (coefficients). From the control surface 
deflection and geometry, the lift coefficient, maximum lift 
angle, moment coefficient, and drag coefficient increments 
are evaluated. The wing lift is defined in terms of lift 
curve slope and maximum lift coefficient. The three- 
dimensional lift curve slope is input directly or calculated 
from the two-dimensional lift curve slope and the wing 
aspect ratio. 

The drag area or drag coefficient is defined for forward 
flight and vertical flight. For small angle-of-attack, the 
drag increases proportional to (1 + K\ a \ x ) , using input 
factor K and exponent X, plus an additional term 
representing separation drag. The induced drag is obtained 
from the lift coefficient, aspect ratio, and Oswald 
efficiency e\ 

(C L -C L0 ) 2 
Dl JteAR 

Conventionally, the Oswald efficiency e can represent the 
wing parasite drag variation with lift, as well as the 
induced drag (hence the use of Qlo)* The wing-body 
interference is specified as a drag area, or a drag 
coefficient based on the wing area. Then 

D = qSC D = qS(C Dp + C Di + C^^) 
is the drag force. 


Figure 11. Wing geometry (symmetric, only right half- 
wing shown). 


Controls 

The wing control variables are flap, flaperon, aileron , and 
incidence. The flaperon and aileron are the same surface, 
generating symmetric and antisymmetric loads 
respectively, hence with different connections to pilot 
controls. With more than one wing panel, each panel can 
have control variables. 


Interference 

With more than one wing, the interference velocity at 
other wings is proportional to the induced velocity of the 
wing producing the interference, with an input factor K ml . 
The induced velocity is obtained from the induced drag. 
For tandem wings, typically K = 2 for the interference 
of the front wing on the aft wing, and = 0 for the 
interference of the aft wing on the front wing. For biplane 
wings, the mutual interference is typically = 0.7 
(upper on lower, and lower on upper). The induced drag is 
then 


c Di = 


(C l -C lo ? 

JteAR 


+ C; 


, ^ ^int a int 
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where the sum is over all other wings. 

The wing interference at the tail produces an angle-of- 
attack change e = E(C L fC La ) , where E = delda is an 
input factor determined by the aircraft geometry. The 
change in orientation of the wing velocity produces the 
interference velocity v£ t at the tail. 

Wing Extensions 

The wing can have extensions, defined as wing portions of 
span b x at each wing tip. For the tiltrotor configuration in 
particular, the wing weight depends on the distribution of 
wing area outboard (the extension) and inboard of the 
rotor and nacelle location. Wing extensions are defined as 
a set of wing panels at the tip. The extension span and area 
are the sum of the panel quantities. 

EMPENNAGE 

The aircraft can have one or more tail surfaces, or no tail 
surface. Each tail is designated as horizontal or vertical. 
The tail is described by planform area S, span b, chord 
c = S /b, and aspect ratio AR = b 2 IS. The tail volume V 
can be referenced to rotor radius and disk area; to wing 
area and chord for horizontal tails; or to wing area and 
span for vertical tails. The geometry is specified in terms 
of S or V, and b or AR or c. 

The horizontal tail can have a cant angle 0 (positive tilt to 
left). The control variables are elevator and incidence. 

The vertical tail can have a cant angle 0 (positive tilt to 
right). The control variables are rudder and incidence. 

FUEL TANK 

The fuel tank capacity (maximum usable fuel 

weight) is determined from designated sizing missions. 
The maximum mission fuel required, (excluding 

reserves and any fuel in auxiliary tanks), gives 

^fuel-cap = ™x(/ fuc , -cap^fuel -miss 5 ^fuel -miss + ^reserve) 

where f^x-cvp a 1 is an input factor. Alternatively, the fuel 
tank capacity W M . cap can be input. 

Auxiliary fuel tanks are defined in one or more sizes. The 
capacity of each auxiliary fuel tank, W aux _ cap , is an input 
parameter. The number of auxiliary fuel tanks on the 
aircraft, N auxtank for each size, can be specified for the 
flight condition or mission segment. Alternatively (if the 
mission is not used to size the fuel tank), the number of 
auxiliary fuel tanks at the start of the mission can be 
determined from the mission fuel. The weight and drag of 
A auxt ank tanks are included in the performance calculation. 


The weight of the auxiliary fuel tanks is an input fraction 
of the tank capacity. The drag area for each auxiliary tank 
is specified, (D/q) ammk . 

PROPULSION 

The propulsion group is a set of components and engine 
groups connected by a drive system. The engine model 
describes a particular engine, which is used in one or more 
engine groups. The components (rotors) define the power 
required. The engine groups define the power available. 
Figure 12 illustrates the power flow. 


POWER POWER 

AVAILABLE REQUIRED 



Figure 12. Power flow. 


Drive System 

The drive system defines gear ratios for all the 
components that it connects. The gear ratio is the ratio of 
the component rotational speed to that of the primary 
rotor. There is one primary rotor per propulsion group (for 
which the reference tip speed is specified); other 
components are dependent (for which a gear ratio is 
specified). There can be more than one drive system state, 
in order to model a multiple- speed or variable- speed 
transmission. Each drive system state corresponds to a set 
of gear ratios. 

For the primary rotor, a reference tip speed V tip _ ref is 
defined for each drive system state. By convention, the 
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“hover tip speed” refers to the reference tip speed for drive 
state #1. If the sizing task changes the hover tip speed, 
then the ratios of the reference tip speeds at different 
engine states are kept constant. 

For dependent rotors, either the gear ratio is specified (for 
each drive system state), or a tip speed is specified and the 
gear ratio calculated ( r = Q dep /£2 prim , Q = V tip _ rcf IR). 

For the engine group, either the gear ratio is specified (for 
each drive system state), or the gear ratio calculated from 
the specification engine turbine speed 
£2 sp^ = (2 jt/ 60) A^ spec and the reference tip speed of the 
primary rotor ( r = Q spec /Q^, = V Up . ref //?). The 

latter option means the specification engine turbine speed 
A/gpec corresponds to V tip _ ref for all drive system states. 

The flight state specifies the tip speed of the primary rotor 
and the drive system state, for each propulsion group. The 
drive system state defines the gear ratio for dependent 
rotors and the engine groups. 

Power Required 

The component power required P com is evaluated for a 
specified flight condition, as the sum of the power 
required by all the components of the propulsion group. 
The total power required for the propulsion group is 
obtained by adding the transmission losses and accessory 
power: P reqPG = P comp + P xmsn + P acc . The power required 
for the propulsion group must be distributed to the engine 
groups. With only one engine group, P reqEG = P reqPG . 

ENGINE GROUP 

The engine group consists of one or more engines of a 
specific type. For each engine type an engine model is 
defined. 

The engine size is described by the power P eng , which is 
the sea-level static power available per engine at a 
specified takeoff rating. The number of engines N eng is 
specified for each engine group. 

Power Available 

Given the flight condition and engine rating, the power 
available P a is calculated (from the specific power SP a 
and mass flow m a ). The flight condition information 
includes the altitude, temperature, flight speed, and 
primary rotor speed; a power fraction f p ; and the states of 
the engine, drive system, and IR suppressor. 

In the engine model, installation losses P loss are subtracted 
from P a ( P av = P a -P loss ), and then the mechanical limit 
applied: P av =mm(P av ,P mechR ). 


The engine model gives the performance of a single 
engine. The power available of the engine group is 
obtained by multiplying the single engine power by the 
number of engines operational (total number of engines 
less inoperable engines): P avEG = (Neng-iVinop)^av The 
propulsion group power available is obtained from the 
sum over the engine groups: P avPG fpPavEG > including 
a specified power fraction f p . 

The drive system rating at the flight condition is 
( Q prim /Q irf) P DSiijmt- Optionally this limit is applied to 
the propulsion group power: P avPG =min(P avPG ,P DS]imit ) . 
Similarly the engine shaft limit at the flight condition is 
optionally applied to the engine group power. 

Performance at Power Required 

The engine performance (mass flow, fuel flow, and gross 
jet thrust) is calculated for a specified power required P q , 
flight condition, and engine rating. The flight condition 
includes the altitude, temperature, flight speed, and 
primary rotor speed; and engine, drive system, and IR 
suppressor states. The engine turbine speed is 

N = r Q 

1 1 eng ' eng aci pnm * 

The engine model deals with a single engine. The power 
required of a single engine is obtained by dividing the 
engine group power by the number of engines operational 
(total number of engines less inoperable engines): 
^req ~ ^ reqEG /(^cng - N in0 p) . In the engine model, 
installation losses P loss are added to P req 

( P — P 4-P \ 

V 1 q 1 req T 1 loss/* 

The performance of the engine group is obtained by 
multiplying the single engine characteristics by the 
number of engines operational (total number of engines 
less inoperable engines): 

^ reqEG ~ (^eng — ^inop ^^req 
^ reqEG ~ (^eng — ^inop ) ^ req ^ ffd 

EG — C^eng — ^inop)^V 

^auxAG = C^eng — ^inop ) ^aux 

The fuel flow has also been multiplied by a factor K^ d 
accounting for deterioration of the engine efficiency. 

Installation 

The difference between installed and uninstalled power is 
the inlet and exhaust losses P loss : P av =P a -P loss and 
P q = P req + P loss . The inlet ram recovery efficiency r] d is 
included in the engine model calculations. The inlet and 
exhaust losses are modeled as fractions of power available 
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or power required. Installation effects on the jet thrust are 
included in the engine model. 

The momentum drag of the auxiliary air flow is a function 
of the mass flow m aux = f am m req : 

^aux — 0 — r hm\ )^aux ^ — (1 — ?7 aux \f ;iux ^‘req V 

where ? 7 aux is the ram recovery efficiency. 

Control and Loads 

The engine orientation is specified by selecting a nominal 
direction in body axes (usually thrust forward), then 
applying a yaw angle, then an incidence or tilt angle. The 
yaw and incidence angles can be connected to the aircraft 
controls. Hence the incidence and yaw angles can be fixed 
orientation or can be control variables. 

The engine group produces a jet thrust F N , acting in the 
direction of the engine; a momentum drag D aux , acting in 
the wind direction; and a nacelle drag D mc , acting in the 
wind direction. 

REFERRED PARAMETER TURBOSHAFT 
ENGINE MODEL (RPTEM) 

Aircraft gas turbine engine performance capabilities are 
formally specified by computer programs known as engine 
decks, which are created by engine manufacturers in an 
industry-standard format. Engine decks are typically based 
on thermodynamic cycle analysis using real engine 
component performance maps. The most important 
performance maps for turboshaft engines are compressor, 
gas generator turbine, and power turbine. These 
component performance maps are critical to obtaining 
realistic off-design engine performance. Design and 
analysis codes calculate aircraft performance for a very 
wide range of operating conditions. This means that 
engine performance must be realistic even far from the 
engine design point. A simple thermodynamic cycle 
analysis that assumes design point component efficiencies 
everywhere is not realistic for such an application. Rather 
than developing models for component performance, a 
model for the total engine performance is used. The engine 
is not being designed. 

The Referred Parameter Turboshaft Engine Model 
(RPTEM) is based on curve-fits of engine performance 
data for existing or projected engines, over a range of 
operating conditions. The curve-fits are typically obtained 
by exercising an engine deck (a computer program). The 
use of referred parameters tends to collapse the data, and 
provides a basis for scaling the engine. The operating 
condition is described by pressure altitude, ambient air 


temperature, flight Mach number, power turbine speed, 
exhaust nozzle area, and either engine rating or engine 
power required. These curve-fits, typically based on real 
engines, are scaled to the required size and adjusted to the 
appropriate technology level to represent a notional 
engine. Engine size is represented by mass flow. Engine 
technology is represented by specific power available and 
specific fuel consumption at maximum continuous power 
(MCP), sea level/standard day (SLS), static (zero airspeed) 
conditions. Engine installation effects (inlet and exhaust 
losses) are also modeled. 

The use of referred parameters to curve-fit engine 
performance data was suggested by David Woodley from 
Boeing during the JVX program (1983). The RPTEM was 
developed and documented by Michael P. Scully and 
Henry Lee of ASRAO, U.S. Army Aeroflightdynamics 
Directorate (AFDD), with a subsequent implementation 
written by Sam Ferguson (1995). 

Operating Environment 

The operating condition and atmosphere give the standard 
conditions (temperature and pressure) for a specified 
pressure altitude; the sea-level standard conditions 
(temperature T 0 and pressure p 0 ); and the operating 
temperature T and pressure p = p sXd . Here the 
temperatures are in deg R or deg K. The engine 
characteristics depend on the temperature ratio 6 = T IT 0 
and pressure ratio 8 = p! p 0 . The flight Mach number M 
is obtained from the aircraft speed. 

The inlet ram air temperature ratio and pressure ratio are 
obtained from M and the inlet ram recovery efficiency 
r) d \ = (1 + 0.2M 2 ) and d M = (1 + r] d 0.2M 2 ) 35 , using 
the ratio of specific heats y = 1.4. 

Engine Ratings 

The power available from a turboshaft engine depends on 
the engine rating. Each engine rating has specific 
operating limitations, most importantly an operating time 
limit intended to avoid damage to the engine. Typical 
engine ratings are MCP (maximum continuous power, no 
time limit), IRP (intermediate rated power, 30 min), MRP 
(maximum rated power, 10 min), and CRP (contingency 
rated power). 

Performance Characteristics 

The engine performance is described by the power 
available P a , at each engine rating and the specification 
engine turbine speed A spec ; the mass flow m and fuel flow 
w required to produce power required P q at engine 
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turbine speed N ; and the gross jet thrust F g at a given 
power required P q . Then the specific power is SP = Pith, 
and the specific fuel consumption is sfc = w IP. 

The reference performance is at sea-level- standard static 
conditions (subscript 0), and MCP (subscript C). 
Referred or corrected engine parameters are used in the 
model: power P/(8 Vfl), mass flow ml(8/'\fd), specific 
power (Plm)l6, fuel flow w/(<5 Vfl), thrust F 18, and 
turbine speed N / V# • 

For each rating R, the performance is characterized by the 
following quantities for sea-level- standard static 
conditions: power Por * specific power SP 0R , and 
mechanical power limit P mQC hR- The mass flow is then 
m 0R = P 0R / SP 0R . The gross jet thrust F g0C is given at 
MCP. These characteristics are at the specification turbine 
speed N spcc . 

The installed power required P req and power available 
P av > P req are measured at the engine output shaft. In 
addition to shaft power, the engine exhaust produces a net 
jet thrust from mass flow that goes through the engine 
core. The fuel flow and mass flow are the total required to 
produce the shaft power and jet thrust. The forces 
produced by mass flow that does not go through the 
engine core (such IR suppressor or cooling air) are treated 
as momentum drag D aux . The relationship between net 
and gross jet thrust is F n =F g -m req V = m req {V j -V) , 
where Vj is the engine jet exhaust velocity. 

The uninstalled power required is P q , the power available 
P a , the gross jet thrust F g , and net jet thrust F n . The 
engine model calculates P a as a function of flight 
condition and engine rating. Installation losses FJ 0SS are 
subtracted from P a ( P av = P a - P loss ), and then the 
mechanical limit applied: P av = min(P av ,P mechR ). 

The engine performance (mass flow, fuel flow, and gross 
jet thrust) is calculated for a specified power required P 
(which might equal the power available), flight condition, 
and engine rating. Installation losses Fj oss are added to 
P ( P = P -i . p } 

1 req V 1 q 1 req 1 loss/* 

Power Turbine Speed 

The shaft power available is a function of the gas power 
available P G and the power turbine efficiency r) t \ 
P a = r] t P G . Generally the power turbine speed N has a 
significant effect on r] t , but almost no effect on P G . The 
model used for the efficiency variation is 

r/ f « 1- 1 (AT/ N opt )~ 1 l* w 


where N opt is the speed for peak efficiency. Alternatively, 
the model parameters are defined at a set of engine speed 
ratios N IN spec . Then the engine performance quantities at 
the required engine speed N are obtained by linear 
interpolation. 

Scaling 

The parameters of the engine model can be defined for a 
specific engine, but as part of the aircraft sizing task the 
parameters must be scaled, in order to define an engine for 
a specified power. In addition, advanced technology must 
be represented in the model. Scaling and advanced 
technology are handled in terms of specific power and 
specific fuel consumption (at SLS static conditions, MCP, 

and W spoc ). 

The engine size is specified as takeoff power P to = P eng : 
power at rating R, for SLS static conditions and 
specification turbine speed N spec . Hence the MCP is 
P oc = P to I r p0R , and the power at all other ratings follows. 

The engine technology parameters SP 0C and sfc oc are 
assumed to vary linearly with mass flow m ^ up to a limit 
m lim , and constant thereafter at SP hm and sfc^. Usually 
the effect of size is that specific power increases and 
specific fuel consumption decreases with mass flow. The 
specific thrust available at MCP is assumed to be constant, 
and the specification power turbine speed decreases with 
the square-root of the mass flow. 

WEIGHTS 

The design gross weight W D is a principal parameter 
defining the aircraft, usually determined by the sizing task 
for the design conditions and missions. The aircraft weight 
statement defines the empty weight, fixed useful load, and 
operating weight for the design configuration. The 
definition of the weight terms is as follows. 

operating weight: W Q =W E +W FUL 

useful load: W UL =W FUL +W pay +W fael 

gross weight: W G =W E +W UL =W Q +W pay +W fucl 

where W E is the weight empty; W F(JL the fixed useful 
load; W pay the payload weight; and W ^ the usable fuel 
weight. The weight empty consists of structure, propulsion 
group, systems and equipment, vibration, and contingency 
weights. If the weight empty is input, then the contingency 
weight is adjusted so W E equals the required value. If the 
design gross weight is input, then the payload or fuel 
weight must be fallout. 
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Table 2. Weight statement (* indicates extension of 
RP8A). 

WEIGHT EMPTY 
STRUCTURE 
wing group 
rotor group 
empennage group 
fuselage group 
alighting gear group 
engine section or nacelle group 
air induction group 
PROPULSION GROUP 
engine system 
propeller/fan installation 
fuel system 
drive system 

SYSTEMS AND EQUIPMENT 
flight controls group 
auxiliary power group 
instruments group 
hydraulic group 
pneumatic group 
electrical group 

avionics group (mission equipment) 
armament group 
furnishings & equipment group 
environmental control group 
anti-icing group 
load & handling group 
VIBRATION (*) 

CONTINGENCY 
FIXED USEFUL LOAD 
crew 

fluids (oil, unusable fuel) (*) 
auxiliary fuel tanks 
other fixed useful load (*) 

kits (*) 

PAYLOAD 

USABLE FUEL 
standard tanks (*) 

auxiliary tanks (*) 


The gross weight W G is specified for each flight condition 
and mission, perhaps in terms of the design gross weight 
W D . For a each flight state, the fixed useful load may be 
different than the design configuration, because of changes 
in auxiliary fuel tank weight or kit weights or increments 
in crew or furnishings weights. Thus the fixed useful load 
weight is calculated for the flight state; and from it the 


useful load weight and operating weight. The gross 
weight, payload weight, and usable fuel weight (in 
standard and auxiliary tanks) completes the weight 
information for the flight state. 

Weight Groups 

Aircraft weight information is stored in a data structure 
that follows SAWE RP8A Group Weight Statement 
format (ref. 18), as outlined in table 2. The asterisks 
designate extensions of RP8A. There are 2 or 3 additional 
levels in the data structure for some weight groups, based 
on the weight breakdown for parametric estimation. 

For each weight group, fixed (input) weights can be 
specified or weight increments added to the results of the 
parametric weight model. The parametric weight model 
includes technology factors. Weights of individual 
elements in a group can be fixed by setting the 
corresponding technology factor to zero. 

The vibration control weight can be input, or specified as a 
fraction of weight empty: W vib = f W foW E . The contingency 
weight can be input, or specified as a fraction of weight 
empty: W cont = f cont W E . However, if the weight empty is 
input, then the contingency weight is adjusted so W E 
equals the required value. 

AFDD Weight Models 

For scaled weights of all components, the rotorcraft 
weight models developed by the U.S. Army 
Aeroflightdynamics Directorate (AFDD) are implemented. 
The weights are estimated from parametric equations 
based on the weights of existing turbine-powered 
helicopters and tiltrotors (and some fixed wing aircraft 
component weights). For some weight groups two models 
are available, designated AFDDnn. Table 3 summarizes 
the statistics of the parametric weight estimation 
equations. 

Figure 13 shows the error of the calculated weight for the 
sum of all parametric weight, accounting on average for 
42% of the empty weight. This sum is composed of the 
structural group (based on the AFDDOO equation for rotor 
blade and hub weights, and the AFDD84 equation for 
body weight), the propulsion group (based on the 
AFDDOO equation for drive system weight), and the flight 
controls group. Based on 42 aircraft, the average error of 
the sum of all parametric weight is 5.3%. The 
corresponding average error is 6.1% for the structural 
group (8.6% for the rotor group alone), 10.9% for the 
propulsion group, and 8.7% for the flight controls group. 
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NDARC SOFTWARE 


Table 3. Statistics of parametric weight equations. 


group 


number of 
aircraft 

average 
error % 

wing 


25 

3.4 

rotor blade 

AFDD82 

37 

7.7 

rotor hub 

AFDD82 

37 

10.2 

rotor blade 

AFDDOO 

51 

7.9 

rotor hub 

AFDDOO 

51 

9.2 

horizontal tail 


13 

22.4 

vertical tail 


12 

23.3 

tail rotor 


19 

16.7 

fuselage 

AFDD82 

30 

8.7 

fuselage 

AFDD84 

35 

6.5 

alighting gear 


28 

8.4 

engine support 


12 

11.0 

engine cowling 


12 

17.9 

air induction 


12 

11.0 

accessory 


16 

11.5 

fuel tank 


15 

4.6 

gear box + 
rotor shaft 

AFDD83 

30 

7.7 

gear box + 
rotor shaft 

AFDDOO 

52 

8.6 

drive shaft 


28 

16.0 

rotor brake 


23 

25.1 

rotary wing 
flight controls 

non-boosted 

20 

10.4 

rotary wing 
flight controls 

boost 

mechanisms 

21 

6.5 

rotary wing 
flight controls 

boosted 

20 

9.7 
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Figure 13. Accuracy of sum of all parametric weight. 


The NDARC program is entirely new software, built on an 
architecture that enables routine extensions and 
modifications. The software has been implemented with 
low-fidelity models, typical of the conceptual design 
environment. Incorporation of higher-fidelity models will 
be possible, as the architecture of the code accommodates 
configuration flexibility and a hierarchy of models. 

The program is written in Fortran 95, using a special- 
purpose software tool to manage the data structures, 
construct the input manual, and automatically generate 
some input and output subroutines. The program has been 
compiled on several platforms and operating systems. 

On typical computers, NDARC execution times range 
from seconds for a job with just a few analysis tasks, to 
minutes for a job that sizes an aircraft based on multiple 
flight conditions and missions. 

Input is in namelist-based text format. The program output 
includes text files formatted for printing and for 
spreadsheets, and special files to support functions such as 
preparing layout drawings. Java graphical user interfaces 
are being developed by the user community to facilitate 
dealing with the input and output. 

The program is supported by complete and thorough 
documentation, including the theory manual (Ref. 6), input 
manual, and data structures manual. A NASA-hosted Wiki 
has been established to support user communication. 

Distribution of the NDARC program is controlled by the 
Software Release Authority at NASA Ames Research 
Center. Source code and documentation are available to 
users, subject to a Software Usage Agreement. 

CONCLUDING REMARKS 

The theoretical basis and architecture of the conceptual 
design tool NDARC (NASA Design and Analysis of 
Rotorcraft) has been described. The principal tasks of 
NDARC are to design a rotorcraft to satisfy specified 
design conditions and missions, and then analyze the 
performance of the aircraft for a set of off-design missions 
and point operating conditions. NDARC provides a 
capability to model general rotorcraft configurations, and 
estimate the performance and attributes of advanced rotor 
concepts. 
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NOMENCLATURE 


Acronyms 

AFDD 

U.S. Army Aeroflightdynamics Directorate 

IRP 

intermediate rated power 

ISA 

International Standard Atmosphere 

MCP 

maximum continuous power 

MRP 

maximum rated power 

RPTEM 

referred parameter turboshaft engine model 

Weights 


design gross weight 


empty weight 

W mo 

maximum takeoff weight 

W SD 

structural design gross weight 

w G 

gross weight, W E +W UL = W Q + W pay 

Wo 

operating weight, W E + W FUL 

W UL 

useful load, W FUL + W pay +W lael 

W m 

payload 


fuel weight 

WfUL 

fixed useful load 

^burn 

mission fuel bum 

W vib 

vibration control weight 

w 

Tr cont 

contingency weight 

Fuel Tanks 

^fuel-cap 

fuel capacity, maximum usable fuel weight 

^auxtank 

number of auxiliary fuel tanks 

w 

rr aux-cap 

auxiliary fuel tank capacity 

Power 

^reqPG 

power required, propulsion group; 
P +P +P 

* comp * xmsn ^ z acc 

^ reqEG 

power required, engine group 

P(ivPG 

power available, propulsion group; 

minCS/pP^XQprim / ^ ref )^D5 limit) 

^avEG 

power available, engine group; 

C^cng — ^inop)^av 
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p 

1 comp 

component power required 

Dlq 

P 

xmsn 

transmission losses 

Aircraft 

P 

acc 

accessory power 

DL 

N- 

mop 

number of inoperative engines, engine group 

^*DS limit 

drive system torque limit (specified as power 
limit at reference rotor speed) 

Aef 

^ES limit 

engine shaft rating 

WL 

^rcf 

Engine 


C AC 

p 

eng 

sea level static power available per engine at 

T 

specified takeoff rating 

c 

N 

’ eng 

number of engines in engine group 

a tilt 


power available, installed; 

Rotor 


— ^loss’ ^mech R ) 

Pa 

power available, uninstalled 

W/A 

Ee q 

power required, installed; P q - P loss 

C w /a 

P, 

power required, uninstalled 

R 

P 

1 loss 

installation losses 

p 

1 mcchtf 

mechanical power limit 

A 

SP 

specific power, Plm (conventional units) 

a 

sfc 

specific fuel consumption, wIP 
(conventional units) 

T 

± design 

m 

mass flow (conventional units) 

r 

w 

fuel flow (conventional units) 


Pn 

net jet thrust 

P 

^aux 

momentum drag 

A 

N 

specification turbine speed 

M at 

V 

Tip Speed and Rotation 

V 

r 

^tip-ref 

reference tip speed, propulsion group 
primary rotor; each drive state 

C T /a 

r 

gearratio; for rotor, 

Q S pec /Q prim for engine 

Pc’ Ps 

^prim 

primary rotor rotational speed, V tip _ ref / R 

^0.75 

^dep 

dependent rotor rotational speed 

0c 0, 

^spec 

specification engine turbine speed 

H,Y,T 

N 

spec 

specification engine turbine speed (rpm) 

M x ,M y 

Axis Systems 

Q 

Pi , P a , f 

I 

inertial 

K 

F 

aircraft 

mean 

Aerodynamics and Loads 


V 

velocity relative air (including interference) 

Wing 

q 

dynamic pressure, %pv 2 

W IS 

a 

angle-of-attack 

s 

P 

sideslip angle 

b 

F 

force 

c 

M 

moment 

AR 

d,y,l 

drag, side, lift forces (aerodynamic axes) 



drag, side, lift force coefficients 



drag area, SC D ( S = reference area) 


disk loading, W D / A ref 
reference rotor area, 2 f A A ; typically 
projected area of lifting rotors 
wing loading, W D IS mf 
reference wing area, 2 S ; sum area all wings 
aircraft control 
control matrix 

component control, c = STc AC + c 0 
tilt control variable 


disk loading, W = f w W D 
design blade loading, W / pAV^a ( V tip = 
hover tip speed) 
blade radius 
disk area 

solidity (ratio blade area to disk area) 
design thrust of antitorque or auxiliary-thrust 
rotor 

direction of rotation ( r = 1 for counter- 
clockwise, r = - 1 for clockwise) 
advance ratio 
inflow ratio 

advancing tip Mach number 
blade flap frequency (per-rev) 
blade Lock number 
thrust coefficient divided by solidity, 
T/pA(QR) 2 a 

longitudinal, lateral flapping (tip-path plane 
tilt relative shaft) 

blade collective pitch angle (at 75% radius) 
lateral, longitudinal blade pitch angle 
drag, side, thrust force on hub (shaft axes) 
roll, pitch moment on hub 
shaft torque 

induced, profile, parasite power 
induced power factor, P t = KP i6Q?a 
profile power mean drag coefficient, 

Cpo = ( a 


wing loading, W = f w W D 

area 

span 

chord, S/b 
aspect ratio, b 2 IS 
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